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Spherically symmetric static empty space solutions are studied in f(R) theories of gravity. We 
reduce the set of modified Einstein's equations to a single equation and show how one can construct 
exact solutions in different f{R) models. In particular, we show that for a large class models, 
including e.g. the f{R) = R — fJ-"^ /R model, the Schwarzschild-de Sitter metric is an exact solution 
of the field equations. The significance of these solutions is discussed in light of solar system 
constraints on f{R) theories of gravity. 



o 
o 

(N 
o 

O; 

in ■ 
(N ; 

> ■ 
cn ■ 

m 

\o ■ 

O ' 
o ■ 

Oh; 

6 



I. INTRODUCTION 

The accelerating expansion of the universe has trans- 
formed our view of the universe from a matter filled cos- 
mos to one dominated by dark energy. Modern day cos- 
mological observations are in contradiction with a matter 
dominated critical density universe whose expansion is 
decelerating and instead we find evidence for expansion 
that is accelerating. Direct evidence supporting the cos- 
mic acceleration comes from the supernovae observations 
p]| and other observations, such as the cosmic microwave 
background and large scale structure 3] , provide more 
indirect evidence. Combining all of the observations, a 
cosmological concordance model has emerged: a critical 
density universe dominated by cold dark matter and cos- 
mological constant -like dark energy. 

The most commonly considered candidate for dark en- 
ergy is the cosmological constant (for a review see e.g. 
Q), but numerous alternative mechanisms for generat- 
ing the cosmic acceleration have been considered. Very 
roughly, one can divide the different alternative explana- 
tions of cosmic acceleration into two classes: those that 
include cosmic fluids with exotic equations of state and 
those that modify gravity. In terms of the Friedmann 
equation one can, again very roughly, consider the for- 
mer to modify the right hand side of the equation, the 
stress-energy tensor, and the latter the left hand side, the 
Einstein tensor. 

Modifications of general relativity (GR) as a source 
of cosmic acceleration have been recently considered in 
numerous works. One particular class of models that 
has drawn a significant amount attention are the f{R) 
gravity models (see e.g. 0, B SH H Ell^ and ref- 
erences therein). These models are a particular class of 
higher derivative gravity theories that include higher or- 
der curvature invariants as functions of the Ricci scalar. 
Such theories avoid the Ostrogradski's instability 
that can otherwise prove to be p roblematic for general 
higher derivative theories 0, h4| . 

A number of challenges have been identified in building 



phenomenologically viable models of f{R) gravity theo- 
ries. Such possible obstacles include instabilities within 
matter outside matter stability of the vacuum 
|18| and constraints arising from known properties of 
gravity in our solar system (see e.g. 0, 0| and 
references therein). In addition, identifying the specific 
functional form of f{R) from cosmological observations 
is problematic since the background expansion does not 
determine /(i?) uniquely [2^ . 

In a number of works, the solar system constraints 
on /(i?) theories of gravity are derived by first confor- 
mally transforming the theory to a scalar-tensor theory 
and then considering the Parameterized Post-Newtonian 
(PPN) limit 24] . This procedure does not seem to 
be without controversy, however H^H^. In this light, it 
is interesting to consider solutions of the modified Ein- 
stein's equations of f{R) theory. Armed with the metric, 
one can hope to study orbital motion directly without 
resorting to conformal transformations. As a first step in 
this direction, we consider vacuum solutions of the mod- 
ified Einstein's equations in this letter. We show how 
one can reduce the set of equations into a single equation 
that one can then utilize to construct explicit solutions. 
As an example we show that a large class of f{R) models 
have the Schwarzschild-de Sitter metric as an exact solu- 
tion (see also for a related discussion). In addition, 
we construct other solutions corresponding to different 
metrics. 



II. f{R) GRAVITY FORMALISM 



The action for f{R) gravity is (see e.g. j27|) 



5* = 



(1) 



where we have set SttG = 1. The field equations resulting 
from this action in the metric approach, i.e. assuming 
that the connection is the Levi-Civita connection and 
variating with respect to the metric g^j/, are 
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where the stress-energy tensor of the gravitational fluid 
is 



1 rl 



fiiy 



with F{R) EE df{R)/dR. 

The standard minimally coupled stress-energy tensor 
T^, derived from the matter Lagrangian Cm in the ac- 
tion (QJ , is related to by 



T;z = r;i/F{R). 



(4) 



In empty space (vacuum), the equations of motion reduce 
to 

F{R)R^, - - V,,V,F{R) + g^,UF{R) = 0. 

(5) 

Contracting these vacuum equations we obtain simply 

F{R)R-2f{R) + 3aF{R) =0. (6) 

This equation is useful, because it allows us to express 
f{R) in terms of its derivatives. If TJJl ^ there is an 
additional trace of stress-energy tensor T™^ in the r.h.s. 
ofEq. ©. 



III. SPHERICALLY SYMMETRIC VACUUM 
SOLUTIONS 

We are interested in spherically symmetric, time in- 
dependent solutions of the empty space field equations. 
From properties of maximally symmetric sub-spaces we 
know that the metric reads as (in spherically symmetric 
coordinates) 
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The (Ol)-component of the Einstein equations is satisfied 
if both p = and s{t,r) is separable with respect to 
its variables. This means that R — and, as in the 
well known case of the Schwarzschild metric, the time- 
dependence can be totally removed from the metric by 
redefinition of time. Here we make the same simplifying 
assumption and thus we consider henceforth only time 
independent solutions i.e. s — s(r) and p — p{r). 



The corresponding scalar curvature is 



R^ —[l-p+{ —){r - - — +77 — 

r^p\ \s p'^ As' 2 s 



(8) 



where we have defined ' = d/ dr. 

Using the contracted equation, Eq. the modified 
Einstein's equations become 

FR^,,-y^V,F ^-^g^,{FR-UF). (9) 

Since the metric only depends on r, one can view Eq. 
© as a set of differential equations for F(r), s(r) and 
p{r). In this case both sides are diagonal and hence we 
have four equations. In addition, we have a consistency 
relation for F{r), 



RF' - R'F + 3(nF)' = 0, 



(10) 



which arises by differentiating the contracted equation, 
Eq. lO with respect to r. Any solution of Eq. must 
satisfy this relation in order to be also a solution of the 
original modified Einstein's equations, Eq. 

From Eq. it is obvious that the combination = 
{F R^^ — V^V^F)/sr^^ (with fixed indices) is indepen- 
dent of the index fi and therefore — A^, = for all 
/i, v. This allows us to write two equations: 

2^+rF'^-2rF" = (11) 

yC A 

F' F' X' 

-is + 4:X - 4rs h 2r'^s' h 2rs 

F F X 

r^s'— + 2r\s" = 0,(12) 
X 

where we have defined X{r) = p{r)s{r). From Eq. I|ll|) 
one can solve for X'/X algebraically and substitute into 
Eq. (inji to obtain X: 



X{r) = s{l + r^- 
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F7) + 



^ 2 // 

— r s ir 
2 ^ 



F 2F + rF' 

F" F' 1 

— ^ —)--r\s". (13) 

2F + rF' F' 2 ^ ' 



Consistency then requires that this form of X{r) satisfies 
Eqs Hll() and (|12(l . giving an equation relating F and s. 
In addition, the modified Einstein's equations give four 
equations relating F and s. However, all of the equations 
have a common factor of the form: 



J 



,(3) 



, -f ArFF' + r'^F'^ - Sr'^FF" 
rF{2F + rF') 
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,8f ^ + r'^F'^ + r^FF'^{4F' + rF") + r^F'^{6F'^ - 3r^F"^ + '(rf (3) - 2F")) + 2rF^{4F' + r{rF^ - F")) 

r'^F^{2F + rF'f 

r^F'^ + r^FF'^(3F' + rF") + PF'^(-3rF"^ + F'(F" + rF^^))) + F^2r(2F" + rF^) - AF') , , 

2s — Tr7:T7Z7. — ^ — ^ = 0- (14) 



r 



^F'^{2F + rF'Y 



Therefore, any pair s(r), F{r) satisfying this equation 
will be a solution of the modified Einstein's equations. 
In addition, if Eq. (|14|l is satisfied, also the consistency 
relation, Eq. pU|l is automatically satisfied. From s and 
F one can then calculate R{r) and in principle construct 
the corresponding f{R) by using Eq. ©. Note that the 
resulting f{R) is not unique due to the presence of an 
integration constant. In addition, a larger degeneracy 
can also exist, e.g. for the SdS-solution discussed below, 
s(r) and F{r) do not determine the f{R) theory uniquely, 
even when discounting the integration constant. 



A. Solutions with constant curvature 

Looking for constant curvature solutions, R = Rq, 
field equations reduce to 



sp + ps 



0, 



1-p 



2p s'^2s ' 2 s 



the 

(15) 
(16) 



which are straightforwardly solvable: 

s{r) 

s(r) = Co H h C2r 

r 



(17) 



where Ci are integration constants. For conventional def- 
initions of space and time we require cq > 0. The scalar 
curvature Eq. © for this solution is i? = 12c2/co. 
Redefining the time coordinate, t i/ \/cq with ci — s- 
ci / ^/cg and C2 — + C2 / , we can always choose cq — 1. 

The Schwarzschild solution in the presence of a cos- 
mological constant, Schwarzschild-de Sitter -spacetime 
(SdS) arising using J{R) = i? + A, has the form 



9ti.y = 



( A{r) 
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-l/A{r) 
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Note that this integral also contains a part of gravita- 
tional energy inside the radius r,„ of the mass distribu- 
tion ^285, while the total energy within the same radius 
also includes vacuum energy, Etot = J d^x (T™° -I- A). 

Comparing the solution, Eq. (|17l) . to the SdS metric, 
Eq. (|18|l , one sees that the two metrics are identical with 
ci = -2Af, C2 = - A/3. From Eq. ® it is clear that this 
metric is a solution for any form of /(i?) for which there 
exists a constant (real) i?o such that i?o/'(^o)~2/(i?o) = 
0. In other words, the SdS metric is an exact solution for 
a set of functions f{R) that satisfy Rof'{Ro)—2f{RQ) = 
such that i?o is real. For example, for the f{R) = R — 
/i^'/i? model, it is easy to see that the SdS metric is a 
solution when R^ = 3/x** = 144c|, i.e. C2 = ^/siJJlM. 
The same exact result also holds for the other commonly 
considered model /(i?) — R — fj,'^/R + eR^. 

The physical interpretation of the parameters M and 
A are not as straightforward as in the case of general 
relativity. The naive identification M = Mc, is prob- 
lematic and a more careful analysis is required. In the 
presence on spherically finite symmetric mass distribu- 
tion the empty space solution we are studying needs to 
be matched to the solution valid inside the mass distri- 
bution at r = r,„. Since the field equations are in gen- 
eral higher order differential equations than in GR, more 
integration constants need to be determined. In partic- 
ular, this means that values of the metric components 
inside the mass distribution depend explicitly on the de- 
tails of the mass distribution, making matching with the 
outside solution non-unique. This can be explicitly seen 
by studying e.g. solutions of spherical shells of different 
thicknesses: the boundary values at Vm depend explicitly 
on the thickness of the shell. This also demonstrates that 
the Birkhoff theorem is no longer valid since the external 
solution depends on the internal mass distribution. This 
is a general property of all (empty space) solutions of 
f{R) theories whenever F differs from a constant. If we 
define the central mass by the gravitational effect it gives 
rise to the external space, the parameter M becomes de- 
fined as the central mass, but it does not coincide with 

dm. 



where A{r) = 1 — 2Af/r— Ar^/3 and A is the cosmological 
constant. The scalar curvature in this case is a constant, 
R = — 4A. For a finite mass distribution, the parameter 
M can in this case be identified as total material mass 



Me 



(19) 



B. General solutions with p{r)s{r) = const. 

From Eq. Hll|) . it is clear that when X = const. = Xq, 
F"{r) = and hence F{r) = Ar + B. Eq. (djl can now 
be solved, giving 



4 



sir) ^ Xo + 



Aci 
2^2 



+ V 



3Br 

SA^B^Xo + 2B^C2 + 2A^ {B^Xq + Aci) In \B/r 
2^4 



(20) 



where Ci are constants. Requiring a SdS-type solution, 
we must choose Xq = 1, which then sets A = and we 
may write also Ci — 6BM, reducing Eq. 1)21) |l to s(r) — 
1 — 2M/r + 02^2 with constant curvature. If we were to 
choose ci = instead, the mass term would be absent. It 
is, however, unclear whether these solutions correspond 
to maximally symmetric (spatial) spaces or to some other 
type of spherically symmetric (but non-singular) cases. 

Requiring that the SdS-type metric is a solution is 
hence equivalent to requiring constant scalar curvature 
and the conclusions of the previous section apply. 

As in the constant curvature case, time can always be 
rescaled so that Xq = 1 . (Alternatively we can choose the 
time scaling of s(r) so that Xq+Aci/2B'^ = 1, which gen- 
erally leads to Xq ^ 1. For examples, see section IlII D\ ) 
Taking Xq = 1 we find that for small values of the radial 
coordinate r, the leading terms of the general solution 
(|20|l read s(r) ^ 1 + Aci/2B^ — ci/3Br where we again 
identify ci = 6BM, leading to the correct form of New- 
tonian potential. However, additional corrections to the 
geodetic motion appear because in general Aci/2B^ ^ 0, 
giving raise to additional parameter constraints. In the 
large r limit the leading contribution comes from the 
-term. It should be noted, however, that s{r) may have 
large finite zeros like in the SdS solution, making the 
limit r —^ oo physically uninteresting. 



C. Asymptotic solutions 

In order to have a better handle on the question of 
uniqueness of the solutions, we consider asymptotic so- 
lutions for which s(r) ^ 1 at large r, mimicking the 
standard Schwarzschild-solution. Our Ansatze are 

s{r) = l-2M/r 
F{r) = For". 

Inserting these into the modified Einstein's equations, 
by requiring that the highest order term vanishes one 
finds that n = or F{r) = const, is the leading term 
in F. Hence there are no new solutions that tend to 
constant scalar curvature in the large r limit along with 
s(r), suggesting that any new solutions will be radically 
different from the Schwarzschild(-de Sitter) solution. 
D. Exact solutions 

We have also found a number of exact solutions to Eq. 
(|14|l . and by considering different Ansatze for s{r) or F{r) 
one can easily find and construct more solutions. A num- 
ber of interesting solutions along with the corresponding 
forms of f{R) are: 



I 

II 
III 



1 - 



2M 



1A.2 



1 - Ur^ 



sor 



Fjr) 
1 - -^r 
For 

For", ?i = 2^i^^^ 



Xjr) 
2 
2 



O Tt*-ra^-4Tt-2 „ „n 



fiR) 



R± 3|jV-i?- 2A + A 



±2y/-R-2A 



2Fo ( 3n(2-n) 
2-n \ rfi-2n-2 



n/2 



R 

-^-2A 

-7J-2A 

3(2-n)Ti 2 
n2-2n-2 



(21) 



These solutions, in particular / and //, could be con- 
sidered as suitable asymptotic limit, either i? ^ or 
i? ^ — CX3, of a more general f{R) having linear f{R) cx R 
term. 



IV. DISCUSSION AND CONCLUSIONS 

We have seen that the set of Einstein's equations re- 
duce to a single non-linear differential equation relating 
s(r) and F{r). In spite of the complicated form, a number 



of solutions is straightforwardly found. The applicabil- 
ity of the general solutions could probably be tested e.g. 
by exploiting parameter constraints appearing from so- 
lar system and comparing these with those arising from 
cosmology. 

Considering the SdS-solution that is present in a large 
class of models, e.g. the R— fi'^/R model with R^ = S/x^ 
or A^ = 3/i*/16, the parameter A can be constrained by 
a number of different observations in the solar system. 
Such are the gravitational redshift measurements, gravi- 
tational time delay measurements by the Cassini space- 
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craft and the perihelion shift of Mercury (see e.g. [2^). 
The tightest constraint arises from the perihehon shift of 
Mercury, for which it is found jl^ that |A| < lO^^^m^'^ 
(the cosmologicahy observed value is roughly lO^^^m"^). 
The solar system observations are hence not able to ef- 
fectively constrain such a metric compared to the cosmo- 
logicahy relevant values. For the other solutions we have 
found, solar system observations are likely to be more 
efficient {e.g. for X 7^ 1). 

In terms of the equivalent scalar-tensor theory, the 
SdS-solution corresponds to a constant field solution. 
It is straightforward to see that the effective scalar 
mass, = V"{(j)), is positive when f'{Ro)/f"{Ro) — 
2f{Ro)/f'{Ro) > 0. This is equivalent to requiring that 
the vacuum state is stable with respect to small pertur- 
bations However, it is not at all clear what, if any, 
role the effective scalar plays since the metric solution 
is independent of the scalar mass. This question will be 
addressed in further work. 

Important questions in addressing the validity of the 
metric solutions presented here are the stability and 
uniqueness of the solution. For example, in order for the 
SdS metric to be physically relevant, it must be stable 
with respect to small perturbations i.e. instead of a test 
mass, one needs to consider the effect of a massive body 
on the metric. Uniqueness of f{R) is also an interesting 
question assessing the physical relevance of a solution. 
Asymptotic considerations indicate that Schwarzschild- 
type metrics lead to constant curvature suggesting that 
new solutions will deviate strongly from the standard 



Schwarzschild metric. 

By considering the boundary conditions of the general 
SdS-solution, we have noted that the Birkhoff theorem 
is not valid for non- trivial f{R) theories as the solutions 
around spherically symmetric mass distribution depends 
on the shape of the distribution. Hence the straightfor- 
ward Schwarzschildian relation between the gravitational 
effect and total energy of a finite mass distribution has 
been broken. By requiring that it holds, we are restricted 
to some class of special mass distributions. Also we are 
lead to ask which distributions are physically relevant 
and what kind of distributions are likely to form from 
collapsing matter. These most interesting but also tech- 
nically extremely tricky questions certainly require fur- 
ther examination because they may offer additional con- 
straints on allowed f{R) models. 

Our results show that in addition to the SdS metric, 
f{R) theories typically also have new different solutions. 
Although further work is needed to determine their phys- 
ical relevance, they offer an interesting new avenue of re- 
search that can guide us in assessing the significance of 
f{R) theories of gravity as a possible solution to the dark 
energy problem. 
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